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Abstract 

Using the recently published model [1] for the collisional energy loss of heavy quarks (Q) in a 
Quark Gluon Plasma (QGP), based on perturbative QCD (pQCD) with a running coupling con- 
stant, we study the interaction between heavy quarks and plasma particles in detail. We discuss 
correlations between the simultaneously produced c and c quarks, study how central collisions can 
be experimentally selected, predict observable correlations and extend our model to the energy 
domain of the large hadron collider (LHC). We finally compare the predictions of our model with 
that of other approaches like AdS/CFT. 

PACS numbers: 12.38Mh 
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I. INTRODUCTION 



High transverse momentum (pr) single non-photonic electrons which have been mea- 
sured in the RHIC heavy ion experiments [2, 3] come dominantly from heavy me- 
son decay The weighted ratio of their p? spectra in pp and AA collisions, Raa — 
daAA/{N c dp^)/(da pp /dp^), where N c is the number of initial binary collisions, reveals the 
energy loss of heavy quarks in the environment created by AA collisions. Initially the az- 
imuthal distribution da/dcj) oc (1 + 2vi ■ cos(0) + 2t> 2 • cos(20)) of light quarks and gluons 
is isotropic and the anisotropy develops during the expansion as an image of the initial ec- 
centricity in coordinate space. The heavy quarks - created in a hard process - are initially 
isotropically distributed in the transverse momentum space. The final v 2 of heavy quarks 
shows therefore how the anisotropy of the light quarks and gluons is transferred to heavy 
quarks and hence reflects this interaction at later times. 

Recently we have published an approach [1] in which we have shown that the energy loss 
as well as the V2{pt) distribution of the single non-photonic electrons in heavy ion reactions 
can be understood in a pQCD based model in which the heavy quarks interact with the 
expanding quark gluon plasma (QGP). In contradistinction to former approaches this model 
has two improvements, 1) a running coupling constant and 2) an infrared regulator in the 
t-channel, which has been determined by physical requirements. 

It is the purpose of this article to explore the details of the interaction between the heavy 
quarks and the QGP in this model, to determine the consequences for observables, to explore 
whether there is a simple way to describe the energy loss in this complicated environment and 
to predict correlations between the simultaneously produced c and c quarks. Furthermore 
we extend the model to the future LHC collider energies and confront the results with other 
theories like the AdS/CFT approach. 

II. THE MODEL 

The model [1] to describe the momentum distribution of heavy quarks or heavy mesons 
produced in ultrarelativistic heavy ion collisions has five major parts which we will describe 
one after the other: 1) the initial distributions of the heavy quarks, 2) the description of 
the expanding quark gluon plasma, 3) the elementary interaction between the heavy quarks 
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and light quarks or gluons, 4) the interaction of the heavy quarks with the plasma and 5) 
the hadronization of heavy quarks into open charm and open beauty mesons. 

A. Initial distribution of the heavy quarks 

For the momentum space distribution as well as for the relative contribution of charm 
and bottom quarks in pp collisions we use the pQCD results in fixed order + next to leading 
logarithm (FONLL) of Cacciari et al. [4]. 
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FIG. 1: (Color online) Transverse momentum distribution of c- and b-quarks in fixed order + next 
to leading logarithm (FONLL) for RHIC ([5]) and LHC ([6]). 

For RHIC energies these results have been published in [5]. They predict a ratio of 
a bbl a cc = 7 10~ 3 . Nevertheless, above pr > Pt cross ~ 4GeV more electrons are produced 
by B-meson decay than by D-meson decay. The uncertainty of p? cross is, however, consid- 
erable. This is due to the uncertainty of the quark masses and of the factorization and 
renormalization scales, fip and fiR. In this work we have taken m c =1.5 GeV, mj, = 5.1 GeV 
and retained the values of and /ip that corresponds to the top curves in the uncertainty 
bands [5] and that are shown in Fig. 1. They provide, after fragmentation into D and B 
mesons and subsequent semi-leptonic decay, the closest agreement with RHIC non-photonic 
single electron data. For LHC energies the initial spectra is stiffer [6] and as well shown in 
Fig.l. The heavy quarks are isotropically distributed in azimuthal direction and therefore 
their u 2 is initially zero. Any observed anisotropy of heavy mesons is due to the interaction 
of their constituents with the medium, and can therefore be used to reveal the strength of 
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this interaction. 

In the E866 experiment at Fermi Lab [7] it has been observed that in pA collisions J / ip 
mesons have a larger average transverse momentum as compared to pp collisions. This 
effect, called Cronin effect, can be parameterized as an increase of < p\ > by 5 ~ 0.2 GeV 2 
per collision of the incident nucleon with one of the target nucleons. We have the option to 
include this effect by convoluting the initial transverse-momentum distribution of the heavy 
quark [5] with a gaussian distribution of r.m.s \Jn co ii(r±) 5o. In this parametrization n co u is 
taken as the mean number of soft collisions which the incoming nucleons have suffered prior 
to the formation of the QQ pair at transverse position f±. It turned out that the Cronin 
effect influences the pt spectra below p? m 5 GeV but is without importance for higher p?. 

In coordinate space the initial distribution of the heavy quarks is given by a Glauber 
calculation. 

B. The expanding plasma 

The expanding plasma is described by a hydrodynamical approach neglecting an even- 
tually existing hard component created by jets. We use the boost invariant model of Heinz 
and Kolb which has been described in detail in [8]. This model reproduces a variety of ex- 
perimental findings. Corresponding to two different equations of state this approach allows 
to calculate two distinct scenarios of the expansion. Either the transition from the QGP to 
the hadron phase is sudden or the system traverses a mixed phase. Hadronization after the 
mixed phase reproduces the spectra of light mesons and is therefore favored by experimental 
data. Without a mixed phase also for heavy quarks the interaction time is too short [1] in 
order to reproduce the energy loss and the azimuthal anisotropy seen in the experimental 
RHIC data. 

Therefore, we use here the model in the mixed phase scenario. We parameterize the 
temperature T(r,t) and the mean velocity u(r,t) of this calculation. This quantities serve 
then to calculate the interaction of the heavy quarks with the medium. They allow to 
calculate the number density of the plasma particles (and hence of the collision rate) as well 
as their momentum distribution. 

At RHIC the initial entropy density for the hydrodynamical calculations is chosen in 
that way that the experimental multiplicity dN c h/dy(y = 0) is reproduced [8]. For the LHC 
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prediction we assume that the soft (thermalized) component contains 1600 < dN^/dyiy = 
0) < 2200. [19] 



C. Elementary interaction between the heavy quarks and the plasma particles 

Using a fixed coupling constant and the Debye mass (mp k gsT) as infrared regulator 
pQCD calculations are not able to reproduce the data, neither the energy loss nor the 
azimuthal distribution characterized by (i^). The novelty of the approach of ref [1] is a new 
description of the interaction between the heavy quarks and the plasma. As compared to 
former pQCD calculations we have introduced a) An effective running coupling constant, 
Oi e s{Q 2 ), determined from electron positron annihilation and non leptonic decay of r leptons. 
b) An infrared regulator in the t channel which is adjusted to give the same energy loss as 
calculated in a hard thermal loop approach. In standard pQCD calculations [9, 10] the gluon 
propagator in the t-channel Born matrix element has to be IR regulated by a screening mass 

a a 

Frequently the IR regulator is taken as proportional to the square of the Debye mass, mo, 
[11] 

H 2 = ml = ^(l + \ n f ) 4vr a s T 2 « (g s T) 2 (2) 
where nj (N c ) are the number of flavors (colors), g 2 s = Anas an d T being the temperature. 

2 

Other approaches use the square of the thermal gluon mass, [12]. In short, jj 2 is not well 
determined. Braaten and Thoma [13] have shown for QED that in a medium with finite 
temperature the Born approximation is not appropriate for low momentum transfer \t\ but 
has to be replaced by a hard thermal loop calculation. Extending their work to QCD we 
have shown [1] that the energy loss, calculated with pQCD matrix elements of the form of 
eq. 1 , only agrees with that calculated in a hard thermal loop approach if /i 2 is much smaller 
and around 

/i 2 « 0.2g 2 s T 2 . (3) 

Employing a running coupling constant and replacing the Debye mass by an effective IR 
cut off (eq. 3) we find a substantial increase of the collisional energy loss which brings for 
the RHIC experiments ^(Pt) as well as Raa{pt) to values close to the experimental ones 
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without excluding a contribution from radiative energy loss. More precisely, the collisional 
cross section has to be multiplied by a K-factor of around 2 (which is assumed to be identical 
for c- and b-quarks) in order to reproduce the data. Thus the difference to the data is of the 
order that we expect for the contribution from radiation energy loss which is not included 
here. 

In this article, we follow the labeling established in ref. [1]. The approach with a running 
coupling constant is dubbed "model E" . In order to point out the influence of the running 
coupling constant we present also some results for the so-called "model C", in which the 
coupling constant is taken as a s (2nT) and fi 2 = 0.15m!,. This model requires K m 5 to 
reproduce the RHIC data. In all calculations, presented here, the corresponding K-factors 
have been applied. 

D. Interaction of the heavy quark with the expanding plasma 

The heavy quarks can scatter elastically with the gluons and light quarks which are 
present in the QGP. The temperature field, determined by the hydrodynamical calculations, 
allows to calculate the density and - together with the local expansion velocity of the plasma 
- the momentum distribution of the light quarks and gluons which scatter with the heavy 
quark. The interaction is described by a Boltzmann equation which is solved by the test 
particle method, applying Monte Carlo techniques. For the collisions between the heavy 
quarks and the plasma particles we apply the elementary pQCD cross sections. We follow 
the trajectory of the individual heavy quarks from creation until hadronization but do not 
pursue that of the plasma particles. The hadronization happens when the energy density of 
the fluid cell falls under a critical value of the energy density. This is 1.64 GeV/fm 3 in the 
scenario without mixed phase and 0.5 GeV/fm 3 at the end of the mixed phase. It is assumed 
that after hadronization heavy mesons do not interact with the hadronic environment. 

E. Hadronization 

The heavy quarks form hadrons either by coalescence or by fragmentation. In our calcu- 
lation the relative fraction depends on pg, on the fluid cell velocity and on the orientation of 
the hadronization hypersurface E as explained below. The coalescence mechanism is based 
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on the model of Dover [14]. To describe the creation of a heavy meson by coalescence we 
start from 



Nq, =D:B = J p Q - p q ■ da : 



d 3 p Q d 3 p q 

2" 



{2nhfE Q (2nhfE q 
xfQ(xQ,PQ)fq(x q ,p q )U(xQ,x q ;p Q ,p q ) (4) 

where Jq et f q are normalized to the number of quarks which go through the hypersurface: 

Pq ■ dax x f Q (x Q ,p Q ) ^ % y E =N Q = 1 (5) 

for hadronization of a given heavy quark and 

J p q ■ da 2 x f q (x q ,p q ) ( 2 ^3fl = N i- ( 6 ) 

/$ is the invariant probability density that a heavy quark at the position xq with momentum 
Pq forms a heavy meson $ with a light quark with x q ,p q , which traverses the hypersurface E. 
f q {x q ,p q ) is the distribution of the light quarks which is assumed to be a thermal Boltzmann- 
Jiittner distribution. Assuming that factorizes we use 

h(x Q ,x q ,p Q , Pq ) = exp ( (X ' ? ~ XQ)2 ~^~ 3;Q),MQ)2 ) x F*(p Q ,p q ). (7) 

uq is the 4-velocity of the heavy quark. Thus in the rest system of the heavy quark is a 
Gaussian function of ||af g — Xq\\- Coalescence requires that in coordinate space the position 
of the heavy and of the light quark are very close and we obtain 

N * = I ^ ' T f^Q,P q )(^Rc) 3 F4PQ, Pq ) , (8) 

J {Znhytjq U Q ■ da 

where da is the unit vector along da. The uq ■ da denominator is positive if the heavy quark 
escapes from the plasma. It counts for the fact that a heavy quark coming out tangentially to 
the critical hypersurface E has a larger chance to encounter its light partner. For F$,(pQ,p q ) 
we take 

\j/ j • (9) 
The coalescence probability is maximal if x q = xq andp q = pq. The normalization condition 

J y^h( x Q, x vPQ,Pq) = !> ( 10 ) 
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where r and p are the relative coordinates between Q and q, relates a\ and R 2 C . We find 

(11) 



AT I \ gggg f rf S f ( n -(^ tt -"^)+^|) « 

N$(xq;pq) = ^ — :/ - u q -da(x Q )e V d / 

Uq ■ da{XQ) Ju q -da(x Q )>0 Uo 

where g q is the degeneracy factor of the light quarks, u q is their 4-velocity and 

5, := hs±<\" x 1 « * (12) 



if ttlq ^> m g . For the calculation we assume a critical temperature of T c = 165 MeV. 
Equation 11 is up to a factor the Cooper- Frye formula which describes the hadronization of 
particles at the surface of the expanding plasma, with an effective inverse temperature /3 e s 
and an effective 4-velocity M ce ii e fr such that j3 e s u ce \\ e s = 3 c u ce \\ + „ UQ a . For a given choice 
of the mass m q , we complete our coalescence algorithm by fixing ad in such a way [20] that 
Nb = 1 for a b-quark at rest in a fluid cell with da = w ce ii, in agreement with the physical 
picture that such a heavy quark can hadronize exclusively by coalescence. The numbers 
and Nb calibrated in this way are thus interpreted as coalescence probabilities, illustrated 
in Fig. 2. 



prob. coal. 




FIG. 2: (Color online) Relative contribution of coalescence of a c(b) quark with a light quark at 
freeze out to the D(B) meson yield as a function of the relative momentum pq/tuq of the heavy 
quark. Heavy mesons which are not produced by coalescence are created by fragmentation as 
described in [5]. 



For momenta above pq = 0.5 GeV the probability to form a heavy meson by coalescence 
falls below one. Because all heavy quarks appear finally as heavy mesons we assume that 
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all heavy quarks which do not coalesce form mesons by fragmentation, as described in [5]. 
As one can see in Fig. 2, high pt heavy mesons are formed exclusively by this mechanism. 

By this hadronization procedure we get a good description of the pt spectrum over the 
whole px range. The physics can best be discussed in terms of Raa which is expected to 
be one if no medium is present. Our results for the D and B mesons for central Au+Au 
collisions at RHIC are displayed in Fig. 3 (for the explication of the different models we 
refer to [1]). In this figure, the upper (lower) limit of the "D-meson" band for model E 
corresponds to m q = 100 MeV (m q = 200 MeV) in equation 11. For B- mesons the difference 
between the two corresponding curves is of the order of the line width. In the following, we 
will retain m q = 100 MeV. 




FIG. 3: (Color online) Raa as a function of pt for D- and B- mesons. We display Raa for two 
parameterizations which describe the experimental data after the results have been multiplied with 
appropriate K- factors (see ref. [1] and section II. C for details). 



III. TOMOGRAPHY AT RHIC ENERGIES 
A. Momentum loss 

We start our analysis with the investigation of the momentum loss the heavy quarks 
suffer while traversing the plasma. In Fig. 4 we display, for central Au+Au collisions at 
yfs = 200 AGeV, the conditional probability density of the transverse momentum loss as a 
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FIG. 4: (Color online) Mean value and variance of the final px momentum distribution of c-quarks 
as a function of their momentum at production for central Au+Au collisions at sfs = 200 GeV. 

function of the initial momentum of the heavy quarks. At high initial momenta (px > 5 GeV) 
we observe a quite broad distribution which narrows for smaller px values. For very low 
initial momenta we see an increase of the transverse momentum during the expansion. If 
their initial px value is smaller than that expected for heavy quarks in equilibrium with their 
environment the interactions with the plasma particles increase their momenta. In Fig. 5 

<p T final [GeV]> 

20 




p T initial [GeV] 



FIG. 5: (Color online) Mean value and variance of the conditional probability density for a c-quark 
with a final transverse momentum px as a function of the initial px- 



the mean value and the variance of this probability density are plotted. Above px = 10 GeV 
we observe, despite of the complex path length distribution, to a very good approximation 
a linear dependence of the px loss on the initial px momentum which can be described 
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by < p|? al >= p™ ltial — 0.08p™ ltial — 5 GeV. Numerically, the constant -5 GeV energy loss 
dominates the —0.08 p™ ltial term even for high p™ ltial but for quantitative comparisons the 
latter is not negligible. This is consistent with the underlying microscopic energy loss, as 
4^ was shown to saturate at large initial momenta due to asymptotic freedom [11]. Also 
the variance depends linearly on the initial px value for high initial momenta. For low 
initial momenta the situation becomes more complex, the final pt approaches there the 
value expected for heavy quarks in equilibrium with their environment. 

To allow for a comparison with other approaches, it is interesting to make the link between 
the energy loss of our model and the transport coefficient 

q = ±- =< kl>ap = « 13 

A v Q At v Q 

which describes the average squared transverse momentum transfer in a single collision 

divided by the mean free A. a, p, B±, At and vq are the heavy quark parton cross section, 

the parton density in the medium, the transverse diffusion coefficient [1], the time between 

two subsequent collisions and the heavy quark velocity, respectively. Fig. 6 shows q as a 

function of the momentum of the heavy quark. 




FIG. 6: (Color online) q as a function of the momentum of the heavy quark for the standard 
parameter set E and for three temperatures of the plasma. [1]. 

B. Dependence of the momentum loss on the creation point in coordinate space 

The momentum loss of a heavy quark depends on the creation point of the heavy quark - 
anti-quark pair. If one wants to know the information about the QGP, which is contained in 
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FIG. 7: (Color online) Left: Average final momentum of the c-quark as a function of its initial 
momentum and of the centrality of its creation point. Right: Relative momentum loss of the 
heavy quark as a function of its initial momentum and of the centrality of its creation point. All 
calculations are done for central Au+Au collisions at = 200 GeV. 

the measured pt spectra, it is important to know from which part of the QGP the observed 
heavy quarks originate. Fig. 7 shows, on the left hand side, the average final pt of the 
heavy (anti) quarks as a function of their initial pt and of the transverse distance of their 
creation point with respect to the center of the reaction, r™. We see three different regimes: 
at low initial pt the average final momentum is independent of the distance to the center. 
These are heavy quarks which come to an equilibrium with the environment. At large 
values of the momentum loss is small because the heavy quarks are too close to the 
(in the hydrodynamical calculation shrinking) surface to interact really with the plasma. 
The third type of heavy quarks are those which have initially a high momentum and have 
been created close to the center. These quarks are really penetrating probes traversing an 
important fraction of the plasma. The momentum loss of these particles is high but it does 
not change substantially between < rip < 4fm. In other words, Fig. 7 shows an explicit 
path-length dependence, despite the rapid decrease of the energy density during plasma 
cooling. This fact contradicts the conclusion of [15]. A complimentary view of the centrality 
dependence of the momentum loss is shown on the right hand side of Fig. 7. There we plot 
the relative momentum loss as a function of the initial transverse momentum p? and of the 
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centrality. The relative momentum loss increases with centrality and with decreasing initial 
momentum. Heavy quarks with a moderate initial momentum suffer the heaviest relative 
momentum loss. There the kinematics of the collisions allows for large angle scattering and 
therefore to a relatively fast approach to equilibrium. The kinematics of the collisions of fast 
heavy quarks with the thermal environment leads only to a moderate momentum transfer 
and therefore the direction of the heavy quark changes only little. Fig. 8 (left) shows the 
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FIG. 8: (Color online) Left: (dN / dpp {&1) ) / (dN / dp^ {itial) ) of c-quarks produced at a [0-2(full), 
2-4(short-dashed), 4-6(dashed-dotted), 6-8(dashed)] [fm] transverse distance to the "center" (sym- 
metry axis) of the reaction for central Au+Au collisions at ^/s = 200 GeV. Right: average trans- 
verse position of the production points of the c-quarks as a function of their final momentum. The 
dashed line corresponds to the px averaged position. The bars mark the statistical uncertainties 
in the simulations. 

dependence of (dN/dp^ 1 ^ /(dN/dp™ ) on the production point of the charmed quark 
- anti-quark pair. Pairs produced at the center of the reaction are highly suppressed at 
large px- Therefore quarks which contribute in this kinematical regime are predominantly 
from the surface and contain little information on plasma properties in the center of the 
reaction. This corona effect can be illustrated alternatively using the correlation between 
the average initial transverse position and the final transverse momentum [21] displayed in 
Fig. 8 (right). There one can see that heavy quarks, passing the hadronization hypersurface 
with px < 3 GeV, originate from larger initial transverse distances than the overall Glauber 
average (~ 3.7 fm). Heavy quarks with a small final transverse momentum come from 
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production points which are more central than the average. We observe a (slow) decrease 
of (r™ ) for large values of pt, because with increasing pf? c the matter becomes less and less 
opaque (see Fig. 7). 

We now discuss a possible way to use the QQ pair as a trigger to probe inner regions of 
the QGP. Within our model, each QQ pair is initially created back to back with the same 
momentum [22]. For the most central production points, the final momentum difference is 
small because the path lengths in the plasma are almost the same for both quarks. The 
more peripheral the pair is produced the more the effective path lengths in the QGP can be 
different. Therefore, the smaller the final pt difference of the simultaneously produced c and 
c pair the more it is probable that it has been produced at a small distance from the center. 
Our approach thus predicts a strong correlation between the final transverse-momentum 
difference of a given QQ quark pair and its initial position in the transverse plane. This 
correlation could possibly be exploited experimentally to discriminate this model against 
other approaches where energy loss is not due multiple independent collisions. In Fig. 9 




FIG. 9: (Color online) Left: Correlation between the average centrality (rip) of the production 
points of the cc pair (labeled iso-contours, in [fm]) and the difference between the final momenta of 
the quarks Ap|? := \pj^ c — PtW f° r various initial momentum pjl = Right: same correlation 
as left for various final pair momenta p^? := T ' c 2 T,c . 

(left), we study further this correlation for simultaneously created cc pairs. We display the 
average transverse distance (r™) of their production points to the center of the reaction as a 
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function of the initial (anti)quark momentum p™, separated for different values of the final 
momentum difference between the c and the c quark. For > 5 GeV, one sees the expected 
correlation: The quarks of pairs created far from the center ((r™) large) have usually quite 
different path lengths and show therefore finally a large momentum difference Apf?. On the 
contrary, for quarks created close to the center the path length is similar and therefore Apf? 
is small. Therefore, by selecting events with small Apf? (w.r.t Pr), one can trigger on more 
central events in order to study their properties. 

In practice, one does of course not have access to p l T as we only measure particles in their 
asymptotic state. Therefore is is useful to study the same correlation as a function of the 
average between p^\ and p^ & i.e. pf, n := T,c 2 T,c . This is done in Fig. 9 (right). Due to 
the energy loss the structure has changed as compared to the left panel. We find now that 
for intermediate values of p^ 1 (g [3,10] GeV), requesting Ap^i 1 ps leads to larger values 
of (r™) , of the order of 5 fm. A refined analysis shows that this small- Apf? crest is due to 
pairs which are produced peripherally and tangentially to the fireball cylinder (p^ _L x"). 
Hence the trajectories of both quarks have approximately the same path-length in matter. 
Although these events are much less frequent than a cc production inside the bulk of the 
QGP, the associated energy loss is rather small, and p^c ~ p™ c , so that they dominate the 
final spectra due to the steeply falling d ^[° d . This interpretation is confirmed by analyzing 

1 dN 1 dN 
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FIG. 10: (Color online) Left: Distribution of the radial distance r™ of the production points of the 
cc pairs for various conditions on their final momenta: no selection (dotted, orange), pj? c > 5 GeV 
(thin, black), pf? > 5 GeV (thick, blue) and p§? c > 5 GeV n Apf? < 0.2 p^ 1 (dashed, red). All 
distributions are normalized to unity. Right: same, with a lower bound of px of 10 GeV instead of 
5 GeV. 

the radial distribution of the creation points, r™, for different conditions on p|? c , as shown 
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in Fig. 10 (left). For large final px values, pfi n c > 5 GeV or > 5 GeV, the corona effect 
is manifest and the central r™ region is clearly depleted w.r.t. the minimum-bias Glauber 
distribution (short dashed). Imposing an additional cut on Apf?, Apf? < 0.2pf, n , we observe 
the disappearance of the corona peak together with a moderate enrichment of the central 
values and an extended "hyper corona" shoulder located at « 6 — 8 fm where those 
tangential emissions take place in which both quarks loose very little energy and therefore 
Ap^ n ~ 0. The right hand side of Fig. 9 shows how one can trigger on more central collisions. 
Increasing values of py increase the sensitivity to central events. With the condition pf, n > 
10 GeV, ApfJ 1 e [1 GeV, 4 GeV], the sample of events is close to that expected from the 
Glauber distribution. The reason for this is that in pQCD calculations with increasing 
energy of the heavy quarks the plasma becomes more transparent, as seen on Fig. 7 (right). 
Energetic heavy quarks produced at small Tt are hence more likely to leave the plasma with 
an appreciable energy and thus they compete in number with quarks produced peripherally. 
Therefore, in the right panel of Fig. 9 < r t > decreases for large pf? 1 . From Fig. 10 (right) we 
see that we nearly recover the Glauber distribution when we apply simultaneously a Apf, n < 
0.2 pf? and ap| n > lOGeV cut. Another consequence of this increasing transparency is seen 

#aa(c) ^aa(c) 




FIG. 11: (Color online) Raa{Pt) ~ here defined has the ratio (dN/dp^J / (dN/dpf) as a function 
of the average of the momenta of the simultaneously produced pair - for different conditions of the 
transverse momentum difference ApJp between the c and the c quark and various selections on the 
transverse distance of the pair creation points with respect to the center of the reaction. The 
left panel shows Raa for all possible Ap^ n and three different bins of r™; the right panel shows the 
result if we apply in addition a cut on the relative transverse momentum of the pair. 

in Fig. 11 and Fig. 12. There, we display Raa(Pt) = (dN/dp^j / {dN/dp l ^j of cc pairs 
for different Apf? selections in comparison with the Raa(pt,c/e) of single (anti) charm quarks 
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(Fig. 12). Only at small pt dN/dfir differs from dN/dp T p due to the Cronin effect, discussed 
in section II. B. Therefore for larger p? Raa{Pt) is a quantity which can be measured. All 
curves in Fig. 11 show a minimum around pt ~ 5 GeV where the relative energy loss 
is most important. The increase beyond px ~ 5 GeV becomes larger if we consider the 
production versus Pt and even larger if we limit the relative transverse momentum of the 
pair. While the second observation is clearly expected, the first one, that Raa behaves 
differently as a function of PT - c + PT > than as a function of pt, c is astonishing and demands 
some explanation in view of the similar profiles observed in Fig. 10. A detailed analysis 
shows that the fluctuations of the average Pt,c ^ Pt ' c are smaller than that of the momentum 
of each of the quarks only. Therefore, the energy loss is less washed out and this explains 
the pt dependence of Raa- 




FIG. 12: (Color online) Raa of c-quarks. Thin lines refers to single quarks, while thick lines 
correspond to the Raa of Pt with and without a Ap T n selection (as in Fig. 11). The full (dashed) 
line refers to model E (C) [1]. 

In conclusion, several relations can be used to validate pQCD based models in general as 
well as our particular model if coincidence data become available because the dependence 
of Raa as a function of Pt for different Ap T n reflects directly the interaction of the quarks 
with the expanding plasma, especially the energy loss as a function of pt- For calculations 
with a running coupling constant (E) the increasing transparency for large Pt quarks is more 
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pronounced than for those with a fixed coupling constant (C). A Apr cut can be used to 
achieve a robust characterization of the energy loss mechanism of heavy quarks. The heavy 
quark pairs are therefore one of the few probes which are sensitive to the expansion of the 
plasma and not only to its properties at the chiral/confinement phase transition. Although 
we have concentrated our analysis on the heavy quarks and not on the observed heavy 
mesons, we expect that the physics seen in Fig. 12 does not change due to hadronization. 
NLO effects at the level of QQ production may modify slightly the conclusions and should 
be included in a future work. 

We see a very complex behavior of the momentum loss of heavy quarks in an expanding 
QGP. It is therefore more than questionable that quantitative predictions of the energy 
loss are possible in models which are based on the average path length of the heavy quark 
trajectory in the plasma. 

C. Azimuthal correlations between simultaneously produced c and c 

Due to the large mass of the heavy quark, interactions between a heavy quark and a 
plasma particle change the direction of the heavy quark only little. We therefore expect 
that 



FIG. 13: (Color online) Distribution of A(f> for different initial pt bins [0-2 (full), 2-5 (dashed), 5-10 
(dashed-dotted), 10-15 (dotted), > 15 (long dashed)] GeV. 

the final azimuthal angle is strongly correlated with the initial one. This is indeed the 
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case for sufficiently high px values, as seen in Fig. 13. There we display the distribution of 
the difference between initial and final azimuthal angle of heavy quarks for different initial 
Pt intervals. The higher px the more small angle scattering dominates and the more we see 
a correlation between initial and final azimuthal angle. There is a sharp transition towards 
an almost flat distribution for p™ ltial < 5 GeV. As already seen, the kinematics allows quarks 
with this initial px to come (almost) to an equilibrium with the environment and therefore 
the correlation is weakened. 



IV. COMPARISON WITH ADS/CFT AT RHIC 

A completely different approach to explain the energy loss of heavy mesons and hence of 
the observed low Raa value at high px has recently been launched by Horowitz and Gyulassy 
[16]. Their model is based on the assumption that QCD is similar to supersymmetric Yang- 
Mills theory and that this theory is dual to string theory in the limit of large N co i or . Whether 
this assumption is justified or not has to be verified. The model allows to calculate the 
momentum loss 

d PT T 2 

— — = —const pt- (14) 

dt M q v 1 

After having implemented this energy loss in a Fokker-Planck approach [18] they could 
predict quite a number of observables which can be confronted with pQCD predictions. One 
of the observables for which the predictions are quite different is the relative energy loss of 
c- and b- quarks. 

pQCD calculations show a much weaker mass dependence (for a given px)- Besides a 
mass dependence of the energy loss in the subdominant u-channel of the gQ — > gQ, which is 

2 

oc the energy loss is only mass dependent for intermediate px {rriQ <C px "C -jr), where 

Q 

4^ oc ln(— ) in the the case of fixed a,. Therefore the difference between the two theories 

ax y rriQ > " 

can be made evident by comparing the Raa (Pt) for bottom and charm quarks. For this 
purpose one may define Rcb{pt) = R c aa(Pt) / Raa(Pt) [18]. In Fig. 14 we compared the 
results of the different theories. It is evident that already the experiments at RHIC energies 
allow to discard one of the theories if high px D and B mesons could be identified. 

Our model yields quite large values of Rcb due to the small value of the IR regulator. 
pQCD calculation with a fixed coupling constant [18] yield smaller values. Due to the 
mild dependence on ttiq, mentioned before, Raa^Pt 1 ) ~ ^AAiPr") as soon as the initial 
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momentum distributions become similar, 4^ ~ t^ 2 -- This is the case for p l £ > 20 GeV/c 

' dpT dpx r 1 — ' 

(see Fig. 1). 




FIG. 14: (Color online) Rcb{pt) = R-aa(pt) / Raa(pt) for different theories. We compare the 
pQCD based "collisional" models C (with K=5, blue) and E (with K=1.8, red) [1] with the 
AdS/CFT calculation for different drag coefficients D/2ttT [16] and for A = 5.5 [16, 17] as well as 
with a pQCD calculation with a fixed coupling constant including radiative collisions [18]. 



V. PREDICTIONS FOR LHC 

Going from the known RHIC to the unknown LHC energy domain we are facing the 
problem that it is not known how the properties of the QGP change with increasing energy. 

We give therefore our results for a range of charged particle multiplicities, 1600 < 
dN c h/dy(y = 0) < 2200, which have been predicted for LHC. We assume furthermore 
that the eccentricity in coordinate space remains the same. Fig. 15 shows the expected 
Raa as a function of pt for D- and B- mesons for model E of [1] which describes best the 
experimental data at RHIC. This model has a K-factor of 1.8. We see that at LHC we will 
experimentally cover the px region in which Raa increases with px- Nevertheless, the Raa 
values are still far away from 1, expected for p T — > oo. 
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FIG. 15: (Color online) Raa for central Pb+Pb collisions at 5.5 ATeV as a function of px for D- 
(left) and B- (right) mesons for model E of [1] for LHC energies as compared to RHIC energies. 
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FIG. 16: (Color online) Rcb{pt) = Raa(pt) / Raa(pt) f° r different theories and for central Pb+Pb 
collisions at 5.5 ATeV. We compare the pQCD based models C (with K=5, blue) and E (with 
K=1.8, red) [1] with the AdS/CFT calculation for different Drag coefficients D/2irT [16] and for 
A = 5.5 [16, 17] and with the pQCD calculation with constant coupling which includes as well 
radiative energy loss [18]. 

The larger pt range at LHC will make it possible to discriminate unambiguously be- 
tween the different energy loss models. Fig. 16 shows Rcb{pt) for three different theories: 
AdS/CFT [16], pQCD with radiative energy loss and constant coupling constant [18] and 
our collisional energy loss model with a K-factor of 1.8 (5) for running (fixed) a s . For mod- 
erate pt (pt < 20 GeV) the spectral form of the c- and b- quarks is different and R~cb{pt) 
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is far from 1, despite the fact that the energy loss becomes more and more similar for c- 
and b- quarks. Above pt = 30 GeV an identical spectral form of the quarks and a con- 
stant energy loss results in Rcb{pt) ~ 1- pQCD calculation are bound to arrive finally at 
values of Rcb{pt) close to one due to the weak mass dependence of the cross section. The 
detailed form of Rcb depends on the cross sections or, more explicitly, on the form of the 
coupling constant and of the IR regulator employed in the pQCD cross section calculation. 
AdS/CFT, on the contrary, predicts even for the highest momenta R cb = 0.2 — 0.3, but it 
has not been demonstrated yet up to which px values the approximations of the approach 
remain valid. 
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FIG. 17: (Color online) V2 for minimum bias reactions Au+Au at 200 AGeV and Pb+Pb at 5.5 
ATeV for model E [1]. The K-factors are noted in the text. Top row: b-quarks (left) and B-mesons 
(right), bottom row: c-quarks (left) and D-mesons (right). 



The azimuthal anisotropy, which has been observed at RHIC energies, will remains visible 
up to LHC energies, as can be inferred from Fig. 17, where we display v 2 for minimum bias 
events separated for b-quarks (top left), B-mesons (top right), c-quarks (bottom left) and 
D-Mesons (bottom right). "e tr m ; n " means hadronization at the end of the mixed phase 
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while "e tr max " means hadronization at the end of the pure QGP phase; in the latter case, a 
K- factor of K = 2-3, has to be applied in order to reproduce the experimental Raa values for 
the most central events. In [1] we found that the experimental values at RHIC can only be 
reproduced when the hadronization takes place at the end of the mixed phase. As it is the 
case for the light hadrons at RHIC, also the i>2 of heavy mesons follows an hydrodynamical 
behavior until pt ~ 2 GeV but the absolute value of V2 is only about half of that of light 
hadrons. 

VI. CONCLUSION 

We have described in detail the predictions of the approach which we have advanced 
recently to describe the energy loss and the azimuthal anisotropy of heavy quarks in the 
environment produced in heavy ion collisions and have extended our calculation toward 
LHC energies. It is based on pQCD calculations with a running coupling constant and 
an IR regulator derived from hard thermal loop calculations. As shown in ref. [1], with 
these new ingredients the energy loss due to elastic collisions is (up to a factor of about 
2) sufficient to produce the observed RaA at RHIC collisions. We have presented several 
observables which allow to test this model. In particular, we predict a large azimuthal 
anisotropy, even at LHC energies and strong correlations between Raa and the transverse- 
momentum difference between the simultaneously produced QQ pair. Correlations between 
simultaneously produced heavy quark pair will allow for triggering on central collisions. 
The identification of D- and B- mesons will allow to reveal whether AdS / CFT describes the 
passage of heavy quarks through matter or whether we are still in the realm of pQCD. 
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